Thermodynamic and diamagnetic properties of weakly doped antiferromagnets 
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Finite-temperature properties of weakly doped antiferromagnets as modeled by the two- 
dimensional t-J model and relevant to underdoped cuprates are investigated by numerical studies of 
small model systems at low doping. Two numerical methods are used: the worldline quantum Monte 
Carlo method with a loop cluster algorithm and the finite-temperature Lanczos method, yielding 
consistent results. Thermodynamic quantities: specific heat, entropy and spin susceptibility reveal a 
sizeable perturbation induced by holes introduced into a magnetic insulator, as well as a pronounced 
temperature dependence. The diamagnetic susceptibility introduced by coupling of the magnetic 
field to the orbital current reveals an anomalous temperature dependence, changing character from 
diamagnetic to paramagnetic at intermediate temperatures. 

PACS numbers: PACS numbers: 71.27.+a, 75.20.-g, 74.72.-h 
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I. INTRODUCTION 

Anomalous normal-state properties of superconducting 
cuprates jl| have stimulated intense theoretical investiga- 
tions of models of strongly correlated electrons describ- 
ing the interplay between antiferromagnetic (AFM) or- 
dering of reference (undoped) insulating substances and 
the itinerant character of charge carriers introduced by 
doping. For the understanding of superconductivity the 
most challenging regime is that of intermediate (opti- 
mum) doping. However, even the apparently simplest 
region of weak doping is not fully understood theoreti- 
cally. 

Recently, the attention in experimental and theoretical 
investigations of cuprates has been given to characteriza- 
tion and understanding of different doping regimes . In 
a simple picture, weak doping should correspond to the 
regime where properties vary linearly with the concen- 
tration of holes, i.e. one can deal with a semiconductor- 
like model where charge carriers (holes) are independent 
and well defined quasiparticles. This requires a nonsin- 
gular variation of thermodynamic quantities with dop- 
ing. However, this scenario has been questioned near 
the metal-insulator transition based also on numerical 
solutions for some model systems Q, e.g. the Hubbard 
model. Alternative possibilities include phase separation 
H , quantum critical behavior M or other instabilities at 
low doping. Still, singular behavior in a planar (2D) sys- 
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tern is expected only at T = 0, while T > should lead 
to a regular variation with doping. 

Among the least understood properties of charge car- 
riers in cuprates and correlated systems in general are 
those related to the coupling of their orbital motion to 
an external magnetic field. Evidently anomalous and not 
understood is the Hall constant in cuprates which re- 
veals unusual temperature and doping dependence [0. 
Another quantity is the diamagnetic (orbital) suscepti- 
bility Xd, which for noninteracting electrons corresponds 
to Landau diamagnetism and seems to be connected 
to the Hall response 0. Anomalous paramagnetic-like 
variation with magnetic field has been noticed within the 
ground state of the t-J model || at low doping. Recent 
T > studies of a single hole within the t-J model || 
confirm the existence of a paramagnetic regime at inter- 
mediate T, though the systems studied were quite small. 
Conclusive experimental results on diamagnetic suscep- 
tibility are lacking ]lQ| , since the orbital part appears 
quite hidden by other contributions, although it could 
be distinguished via the anisotropy. 

The aim of this paper is to study the thermodynamic 
properties and orbital response of correlated electrons at 
finite temperature in the low-doping regime. Most nu- 
merical studies of the t-J model have so far focused on 
the ground-state properties [Q, employing exact diago- 
nalization of small systems, projector Monte Carlo, and 
density matrix renormalization group [Jl2j (DMRG). Re- 
cently, the finite-temperature Lanczos method (FTLM) 
has been introduced, which allows insight into the stat- 
ics and dynamics at T > 0. In previous applications 
certain thermodynamic quantities have also been inves- 
tigated as a function of doping. In this paper we focus 
on the low doping regime, where the method can be com- 



2 



pared with the alternative approach, a novel adaptation 
of the worldline quantum Monte Carlo (QMC) cluster 
method which allows for the study of much larger 
systems at least for temperatures T > T_ below which 
the minus-sign problem sets in. Large systems are partic- 
ularly important for the study of diamagnetic response 
which appears to be quite sensitive to finite size effects. 
In both cases, new ways of dealing with the magnetic 
field are introduced. Related QMC methods have been 
used to study nonmagnetic properties of the t-J model, 
in an exploratory calculation for doped chains and for 
ladders with 1 and 2 holes p5[ , in two dimensions at 
J — > with 1 or 2 holes Eq^ and for chains at finite J in 
a background of no holes ||14| . 

In the following, the planar t-J model as a representa- 
tive model for strongly correlated electrons and electronic 
properties of cuprates is studied, 



H = -t (A ~, 



H.C. 



(ij) 



(1) 



where cj CT , c ia are fermionic operators, projecting out 
sites with double occupancy. To approach the regime of 
strong correlations close to the real situation in cuprates, 
J/t = 0.4 is used in most numerical calculations. We also 
use fcs = h = 1. 

The paper is organized as follows. Section II of the 
paper is devoted to a brief introduction of both numerical 
techniques employed, QMC and FTLM. In Sec. Ill results 
for several thermodynamic properties in the low-doping 
regime are presented and discussed. Sec. IV is devoted to 
the discussion of the orbital susceptibility of the system. 



II. NUMERICAL METHODS 

Results are obtained independently by the worldline 
QMC method and the FTLM. Wherever possible, re- 
sults of both methods for doped systems are compared 
and presented relative to the undoped Heisenberg AFM. 
For large enough systems we expect to reach a typical 
behavior in the low doping regime. 



A. Worldline quantum Monte Carlo method 

The loop cluster algorithm (LCA) for the world-line 
QMC has been introduced by one of the present authors 

"Tel. 



fli3| and recently adapted also to the t-J model |15 



We briefly describe the worldline representation of the 
quantum QMC. The Hamiltonian, Eq. |l|, on a 2D square 
lattice can be split within the standard Trotter-Suzuki 
decomposition |lj], into four parts H = Hi + H 2 + 
H3+H4 consisting of mutually commuting terms. This is 
equivalent to the well known checkerboard decomposition 



of Hamiltonians in ID. The partition function is 

Z = Tre-^= hm t^-P^+hm+h^m ^ 

M-> 00 ' 
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where $ = (3/M and j3 = 1/T. The summation is taken 
over the complete orthonormal set of states \4>i). Within 
each imaginary time step (3 the time evolution operator is 
applied. Since the Hamiltonian is total spin conserving, 
we can track time evolution of a particular spin along its 
so called worldline (WL). Because of the cyclic property 
of the trace, the WLs are periodic in the imaginary time 
interval [0, /?]. The time evolution operator acts only on 
2x2 plaquettes and the weight of the configuration W(C) 
factorizes into a product of plaquette weights. The par- 
tition function 



C C pGC 



(3) 



is formally that of a (2 + l)-dimensional classical system. 
The thermal average of an observable O can be obtained 
by 



<O> = ~XV(C)0(C)- 



(4) 



Such thermal expectation values are calculated by means 
of Monte Carlo (MC) importance sampling, where a 
sequence of configurations Ci (Markov chain) is con- 
structed, which obeys detailed balance and reproduces 
the correct Boltzmann distribution W(C)/Z. Thermal 
expectation values now become simple averages 



(O) 



lim — 
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E 
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(5) 



In practice, Monte Carlo runs are finite, K < 00, lead- 
ing to statistical errors which can be calculated from the 
standard deviation of partial data sets |L3| . 

In standard local algorithms an update from one con- 
figuration C to another C in the Markov chain represents 
a small local change of the WLs. Therefore, consecutive 
configurations are highly correlated, which drastically in- 
creases the necessary number of Monte Carlo steps. Such 
difficulties are overcome in the LCA [[l3] which introduces 
global (nonlocal) stochastic updates that effectively re- 
duce the correlations. In the LCA formulation also the 
continuous time limit (3 — * can be taken J 19 avoid- 
ing the second order systematic error of Eq. |2|. For cer- 
tain observables improved estimators can be easily con- 
structed allowing a potential reduction of statistical er- 
rors. For more details we refer to the introductory paper 

& 
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The LCA has recently been adapted to the t-J model. 
The update procedure is split into three substeps, allow- 
ing the application of the standard LCA for the S =1/2 
antiferromagnetic Heisenberg model or for free fermions 
in all three cases. Within each substep, only updates be- 
tween two of the possible three states (f, |, and hole o) 
are performed. For the weights of particular plaquettes 
and other technical details we refer to [p"5| . 

In case of negative weights W(C) < 0, their magnitude 
|W(C)| is taken for construction of the MC procedure, 
since the negative W{C) cannot be taken as a probability. 
Eq. (Q) becomes 



(0) = 



(sign • O) 



w\ 



(sign) 



W\ 



(6) 



where 



denotes the expectation value with respect 



to the absolute value of the weight. In systems with 
such a "sign problem", the average sign (sign), w , often 
becomes exponentially small with increasing system size 
and decreasing temperature T, leading to a blow up of 
statistical errors Ji"8|. 

Let us briefly comment on the origin of negative signs 
in the WL formulation of the t-J model. In the system 
with no doped holes the only source of negative weights 
are plaquettes, where two opposite spins exchange their 
positions representing a spin flip. Because of the period- 
icity of WLs in time direction and the absence of holes 
spin flips always occur in even numbers, producing no 
net negative sign. In the pure Heisenberg model, this 
sign can also be transformed away by rotation of spins on 
one sublattice, resulting in all-positive plaquette weights 



For one hole doped into the AFM one would naively 
not expect a sign problem, e.g., in this case there is no 
sign in the exact diagonalization approach. Examining 
the particle WLs surrounding the hole WL one finds, 
however, that an exchange of two fermions can occur 
when t ^ and J ^ 0, producing an odd number of spin 
flips, i.e. a negative sign, as can be seen schematically in 
a small 2x2 system, 
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where a loop motion of the hole Q around the system 
and a consecutive spin flip <~~> reproduce the original 
configuration with two fermions f and | interchanged. 
Measuring the sign here reduces to spin flip counting. 

For higher concentration of holes a more general ex- 
pression for the sign of the configuration can be obtained. 
It links fermion WL (perm/) and hole WL permutation 
(pernio) 

sign(C) = (-1)P C ™/ = (_1)P«™* . sign(H/(C)), (7) 
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FIG. 1: Average sign in QMC for different values of the 
anisotropy 7 and number of holes Nh for a 6 x 6 system with 
J = 0.4t 



the use of the improved estimators since for every observ- 
able a separate algorithm must be devised. 

To follow the emergence of the sign problem as well 
as the development of diamagnetic properties it is con- 
venient to generalize the isotropic spin interaction term 
of the model Eq. ([j]) to an anisotropic one with general 
anisotropy parameter 7, 



Hi 



(ij) 



SJS, 
3 



(8) 



This modifies the pure spin substep (|, J.) of the t-J LCA 
so that otherwise independent loops are "frozen" together 
Jl3] | into clusters and updated stohasticaly. 

The results for (sign) as a function of inverse temper- 
ature (3t = 1/UbT are presented in Fig. [I]. Note that for 
a single hole in the system, Nh = 1, the relevant tem- 
perature scale in the anisotropic case is T_ ~ 7J, i.e. 
there is no sign problem for 7 = 0. At T < T_ the sign 
starts to deteriorate rapidly, as can be seen in Fig. ^, pre- 
venting the investigation of low temperature properties. 
As expected (sign) decreases by adding additional holes 
Nh > 1. For this reason, within the doped t-J model only 
chains and coupled chains |H| have been investigated by 
LCA so far, and in two dimensions the limit J — > with 
Nh = 1,2 |L6]]. Recently, though, a way around the sign 
problem has been proposed by calculating fermion prop- 
agators for a background of no holes H] . 

For a fixed number of holes, the average sign will con- 
verge as the system size increases. This convergence 
could be taken as a another criterion that the limit of 
a dilute system has been reached. 



B. Finite-temperature Lanczos method 



so that for low doping it is preferable to measure perm/j 
rather than permj. The sign problem also complicates 



In the analysis of the t-J model the exact diagonal- 
ization of small systems using the Lanczos algorithm has 
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been extensively employed [ fLl| , predominantly in the in- 
vestigation of the static and dynamic properties of the 
ground state. More recently a FTLM combining the 
Lanczos procedure and random sampling was introduced 
pl| . p2{ , allowing the calculation of T > static and dy- 
namic properties of correlated systems. The application 
is particularly simple for an arbitrary function of con- 
served quantities, e.g. f(H,S z ), 

z « ^EE^iwdi 2 , 0) 

n=l i=0 
K M-l 

(/) » ^fEE/^"^r)e-^"iHv-r>i 2 ,(io) 

n=l i=0 



3 







AE 1 FTLM N=20 






---- AE 2 FTLM N=20 


- 




o AE, QMC iV=6x6 


- 




• AE 2 QMC /V=6x6 


k. 



















0.0 0.5 1.0 1.5 2.0 



where \tp"), E™ are (approximate) eigenfunctions and en- 
ergies, respectively, obtained by diagonalization within 
the reduced orthonormal set, generated from the initial 
functions \n) in M Lanczos steps. N st is the dimension 
of the complete basis. K initial functions |n) are chosen 
at random but with good quantum number S z . Usually 
it is enough to choose M, K <C N st . For a more detailed 
discussion of the method and results we refer to |^] . 

It is expected that T > reduces the finite-size effects 
of the measured quantities. It is however important to re- 
alize that for a particular system, finite-size effects start 
to be pronounced at T < Tf s where, e.g., some charac- 
teristic length-scale becomes larger than the system size. 
In our case of low doping, N < 20 and J = 0.4 t, we 
find Tf s ~ 0.4 J. All our results are presented for T > Tf s 
where Z(T is ) ~ Z* . In the present study Z* = 30 so 
that at least 30 states are sampled in the thermal aver- 
ages Q . It should be stressed that the FTLM gives also 
the correct ground state within the chosen small system. 



III. THERMODYNAMIC PROPERTIES 

Thermodynamic properties (D(ch) depend on the hole 
concentration Ch — Nh/N. For the weak doping limit, 
one would expect a linear dependence for most quantities. 
For a finite system size, the relevant parameter is thus 
the number of holes Nh doped into the AFM. In the low- 
doping regime it makes sense to represent the results as 
a difference 



Ad = 0(N h =i) - 0(N h =i - I). 



(11) 



To distinguish the change in a particular quantity with 
doping, this notation is used in the following. If, e.g., 
A02 behaves quantitatively as AO\ one can conclude 
that the quantity changes linearly with the number of 
added holes Nh, i.e. the holes behave as independent en- 
tities, and the system sizes are large enough so that the 
low-doping regime has indeed been reached. Such be- 
havior is however not the only possibility at low doping, 
since one can expect, e.g., even-odd effects in the case of 
pairing of holes. 



FIG. 2: Internal energy per hole AE vs. T as calculated within 
the FTLM and the QMC method for J = OAt. ■ represents 
the one-hole ground state result for Af = 4x4 from Ref. |nj . 
A and T are our ground state results for AEi and AE2 for 
N = 20, respectively. 



A. Internal energy, specific heat and entropy 



The internal energy, defined by 



E = 



df3F 
~~d/3~ 



1 dZ 
Z~di3' 



(12) 



is calculated within FTLM as (E) in Eq. (|i"o| ) and in 
QMC as an expectation value of the corresponding oper- 
ator, E — Ylc E(C)W(C). Results of both methods are 
presented in Fig. 0. 

From Fig. || we first conclude that for AEi the re- 
sults at TV = 20 (which overlap also with the results at 
N = 18) obtained via the FTLM are essentially equiva- 
lent with QMC results for much larger lattices, at least 
for T > 0.2 t reached by the QMC method. We note also 
a close agreement between QMC AE2 and AEi , confirm- 
ing the assumption of the low-doping regime and holes 
as independent quasiparticles in the T window presented. 
In FTLM results, on the other hand, the difference be- 
tween AE2 and AE\ is already visible, since Nh = 2 
here means already an appreciable doping Ch = 0.1. For 
T = the difference AE2 — AE\ equals to the binding 
energy jllj] and in the continuum corresponds to the sec- 
ond derivative of the ground state energy with respect to 
the doping. In the chosen parameter regime (J = OAt) 
the binding energy is negative and thus pointing to the 
attractive interaction between the holes. With the in- 
crease of the temperature the bound state disintegrates 
and the difference AE2 — AE\ approaches zero. In the 
case of a small system the difference can even become 
positive but vanishes with increasing system size. 

It is also evident from Fig. | that AE(T) is not a 
monotonous function. The ground state of a single 
hole introduced into the AFM is quite well understood 
via analytical approaches |23 and numerical calcula- 
tions Q. For J/t = 0.4, the zero temperature result 
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FIG. 4: Entropy increase AS 1 (in units of fcs) for a single 
hole obtained via the FTLM and the QMC method, on dif- 
ferent size lattices and J — OAt. Long-dashed lines represent 
AS(oo). 



pears in AC as a pronounced dip which slightly weakens 
and shifts its energy scale J to lower values with doping, 
as can be seen from the line-shape in Fig. |^. 

The entropy is 



FIG. 3: Specific heat per site C/N for undoped system (top) 
and change with doping AC (bottom, both in units of fcs) 
vs. T as calculated within the FTLM and the QMC method 
for J = OAt. 



S = (3(E-F) = /3E + lnZ 



We reconstruct it from the specific heat 



(14) 



AE(0) ~ — 1.444 can be explained well by the interplay 
between the gain of the kinetic energy represented by the 
hopping term H t and the loss of local AFM correlation 
energy around the hole. 

AE(T) has not been considered so far. An interpreta- 
tion of its behavior can be given as follows. Introducing 
a single hole into an AFM destroys the local AFM spin 
order and thus increases the exchange energy. The in- 
crease is however expected to disappear at T > J where 
the spin system becomes disordered. On the other hand 
the ground-state kinetic energy in a disordered spin sys- 
tem is quite similar to the one in an AFM, hence the 
decrease of the internal energy AE for T > J. This re- 
mains valid for T < t where also higher hopping-related 
states become populated and finally AE(T — > oo) — > 0, 
explaining turn back up for T > 0.74. 

The specific heat defined by 



C'-^-B 2 
6 " dT ~ * 



1 d 2 Z 

Z d(3 2 



1 dZ 
Z~d{3 



(13) 



is obtained as P 2 [(E 2 ) - (E) 2 } within the FTLM and the 
QMC method. The results are presented in Fig. ||. 

The main effect of introducing holes into the AFM in- 
sulator on C is to decrease the peak at T ~ J. This ap- 



C - 1 dT 



(15) 



by numerical integration from high temperatures T ~ oo 



f T AC 

AS(T) - AS(oc) = / — dT 



(16) 



The high-temperature integration constants are chosen 
so that AS(oc) = AlniV st . 

In contrast to AE and AC discussed previously, AS 
is not linear in in the low doping limit. In analogy 
with low-concentration systems, like the dilute classical 
gas, entropy is expected to scale as AS oc \ A(Nch hic/j)|, 
i.e. A Si oc InN, so the change still depends explicitly 
on the system size. This is also realized in FigfJ, where 
ASi still varies with N, yet curves for different TV appear 
parallel down to the lowest reachable temperatures. 

It is particularly remarkable how large the entropy in- 
crease ASi ^ 1 is even at the lowest T < J. This 
is indeed consistent with AS measured in cuprates Q. 
While this could be attributed partly to the logarithmic 
dependence on Ch, at the same it is evident that the be- 
havior is much closer to a system of classical particles 
than to a degenerate electron gas. 
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FIG. 6: Effective square moment (Curie constant) per hole 
A (S 1 ^) vs. T, obtained from the data in Fig. E| For compari- 
son also the result for noninteracting classical spins (NCS) is 
presented (dashed line). 



at the lowest reachable temperature T ~ J/2 each hole 
effectively adds just one spin, i.e. A^S 2 ') ~ 1/4- 



IV. ORBITAL SUSCEPTIBILITY 



FIG. 5: Spin susceptibility for undoped system y s (top) and 
change of the spin susceptibility with doping Ay s (bottom) 
vs. T obtained within the FTLM and the QMC method for 
J = 0.41 



B. Spin susceptibility 

The uniform spin susceptibility can be evaluated as a 
thermodynamic quantity from 



Xs 



N 



(17) 



where S z = Sf is the conserved total spin. In the 
FTLM then Eq. (10) can be applied, while within the 
QMC method %s is related to the number of spin up and 
down WLs. 

It is instructive to present results both for A\ s with 
respect to the undoped AFM, Fig. ||, as well as for the ef- 
fective Curie constant (difference of the square moment) 
per hole A (S 2 ) = NAy s /P in Fig. | 

The results in Figs, pi ra are easy to interpret for high 
T > t. Each hole introduced into the system reduces 
the effective Curie constant by one spin, i.e., A(S 2 ) = 
— 1/4. On the other hand, at low T < J the situation is 
reversed since A\ s > 0. This increase can be attributed 
to the relaxation of the AFM order by the hole doping. 
Note that in an AFM, Xs achieves a maximum at T ~ J 
while below that temperature it is reduced due to the 
longer range AFM order. It is interesting to note that 



In order to investigate the orbital response of the sys- 
tem, a homogeneous magnetic field B perpendicular to 
the plane has to be introduced. When we discuss the or- 
bital magnetization and susceptibility, B enters only in 
the kinetic term of Eq. (oh, via the Peierls construction 



-*£(■ 



H.c), 



(18) 



where the phases are given within Landau gauge as 



-A.:r,;.-R, 



B(0,x,0), (19) 



with R, 



Yi. The relevant parameter for the 



strength of B is the dimensionless flux per plaquette 
a = 2itBa 2 / <j)Q, where a is the lattice spacing and 
(pa = h/e is the unit quantum flux. 

The dc orbital susceptibility of the system in the ex- 
ternal magnetic field is 



Xd 



-Mo 

Xo 
" (3 



d 2 F ^ e 2 a 4 d 2 F 



dB 2 



1 d 2 Z 



h 2 da 2 



1 dZ 
Z~da 



(20) 



where xo — Moe 2 a 4 / h 2 . So far \d has been investigated 
only for a single hole = 1 by high-temperature ex- 
pansion at J = 0, and by the FTLM for J > (9). In the 
latter study it was realized that results are quite sensitive 
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to finite-size effects, so it is desirable to get correspond- 
ing results also via the QMC method, where much larger 
lattices can be studied. 

Let us first derive the expression for the orbital suscep- 
tibility within the QMC method. As seen from Eq. ( Km , 
the magnetic field affects only the hopping of the elec- 
trons. In the WL representation this concerns matrix 
elements (<f>\ e~@ H \<fi') in Eq. (||). For the plaquette rep- 
resenting the hopping between sites i and j the matrix 
elements become 



-te l > 
-te- t9 *i 



(21) 



written in the |o f), |f o) base. The imaginary time prop- 
agator in the same base is 



ch pt e ie « sh fit 
i6 » sh.pt ch pt 



(22) 



Thus the plaquette weights along the hole WL obtain 
an additional phase factor W(p,a) — W(p)e l8 ^ p \ The 
weight of the whole configuration is a product, Eq. (||), 

W(C, a) = J] W(p, a) = W(C) ]J exp (i6(p)) (23) 
pec pec 

where the phases sum up 



exp 



[i^2 e (p) = ex P if 



(r)dr 



exp 



i — <p A(r)dr 



(24) 



The integral runs along the hole WL. Here S is defined 
as the oriented area of the hole WL projected onto the 
plane in units of the lattice plaquette area a 2 . For more 
holes, S generalizes similarly to the sum of all hole WL 
areas. 

Now we can write the partition function in the mag- 
netic field as 



Z = J2w(C,a)=Y, e 



iaS(C) 



W(C). 



(25) 



For a given configuration C there always exist C (imag- 
inary time inversion) with the same weight but S(C) = 
—S(C), therefore the exponential in Eq. fl25|) can be re- 
placed by a cos function. For B = we have (S) = and 
obtain the zero-field susceptibility from Eq. (pOl), 



Xd 



-Xo 



<£!> 

p 



(26) 



Xd can be thus measured without the presence of a mag- 
netic field. This is just another consequence of the more 
general fluctuation-dissipation theorem. Even though S 2 
is strictly positive, the thermal average in Eq. (|2^) can 
become negative because of correlations between the sign 
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FIG. 7: Finite-size scaling of \d from QMC with one hole, for 
different T and J = 0.4i. Solid lines serve as a guide to the 
eye and are fits to the 1/N 2 and 1/N 4 dependence. 



of the weight and the area S of the hole WL. From Eq. (§) 
we can deduce that (S 2 ) < when the configurations 
with negative sign tend to have larger 5 2 (C)|IF(C)| than 
configurations with positive ones. 

The hole WL can obtain nonzero spatial winding num- 
ber due to the periodic boundary conditions and small 
system size. E.g., the WL can run along the imaginary 
time, cross the system boundary, and complete time pe- 
riodicity reconnecting with its spatially periodic image. 
In that case the area S as defined by Eq. ([m]) has no 
physical meaning. Therefore we restrict our simulation 
so that only the zero spatial winding loop updates are 
generated (a discussion on fixing the winding numbers 
can be found in |2j|). The hole is allowed to cross the 
system boundary as long as it does not increase the wind- 
ing number. The effect of the restriction is analogous to 
the movement of the hole doped into an infinite periodic 
spin background with a unit cell equal to the size of the 
system. The results for thermodynamic quantities pre- 
sented in the previous section agree within the errorbars 
with the unrestricted case. This restriction is weaker 
than closed boundary conditions, resulting in smaller fi- 
nite size effects. 

The introduction of finite B > into the model, 
Eq. (|l8|), reduces the translational symmetry and thus 
for a given system size increases the required minimal 
base set used in FTLM. In the present study a few mo- 
bile holes on a system of tilted squares with N up to 20 
sites and periodic boundary conditions are considered. It 
is nontrivial to incorporate phases 9ij corresponding to 
a homogeneous B, being at the same time compatible 
with periodic boundary conditions ^6|. This is possi- 
ble only for quantized magnetic fields B = mBo, where 
B = <f>o/N. 

Xd from Eq. ( p6[ ) can be calculated in FTLM only by 
taking a numerical derivative of the free energy F = 
— TlnZ, with Z from Eq. (Q). Finite systems provide 
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FIG. 8: Orbital susceptibility \<i vs. T obtained via the QMC 
method and the FTLM method for J = OAt. Regarding the 
FTLM xi2 for N = 20 should be most relevant (thick line) . 



F(a) only for discrete values of a. \d can be obtained 
by fitting the F(a) dependence to the parabolic form 
F(a) = F(0) + co? in two points a = i ■ 2tt/N and 
a = j ■ 2ir/N. The corresponding results for the suscepti- 
bility are denoted by Xij ■ I n small systems (TV < 20) the 
introduction of B > can lift some zero-field degenera- 
cies. The xoi values are therefore systematically affected 
by larger finite-size effects. For the system with N = 20 
both values Xoi and X12 agree quite reasonably with the 
QMC data. 

Let us first comment on the validity of results for Xd- 
Since Xd deals with an orbital current represented by loop 
motion of charge carriers (holes), it is much more sensi- 
tive to finite-size effects || than most other correlation 
functions. This was also the main motivation to employ 
the QMC method, where much larger systems can be 
reached. In Fig. |?], finite-size scaling is performed for the 
QMC data for Xd(T) for the case of Nh — 1. We can see 
that the different T points do not cross upon changing 
the system size N. Thus, at least qualitatively, results 
do not depend on the system size. Therefore, the sizes of 
choice for the QMC systems will be 6 x 6 and 8x8. 

In Fig. H Xd obtained via both methods is presented. 
For T>f, the response is diamagnetic and proportional 
to T~ 3 as well as essentially J-independent ||. The 
most striking effect is that the orbital response below 
some temperature T p turns from diamagnetic to param- 
agnetic, consistent with the preliminary results obtained 
via the FTLM ||. In order to locate the origin of this 
phenomenon, results for different J and anisotropies 7 
are shown in Fig. ^. It appears that T p scales with 7J, 
i.e. at J = the response is clearly diamagnetic at all T, 
and for 7 = 0, J > no crossing is observed with either 
method. 

At lower temperatures T < Td <C T p , the diamagnetic 
behavior is expected to be restored. This follows from the 
argument that at T — > a hole in an AFM should behave 
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FIG. 9: Xd vs. T for N h = I and different J and 7. QMC 
results for the 6x6 system are shown with different dots and 
are labeled with appropriate values of J and 7. Solid lines 
correspond to the N = 20 FTLM for the same values of J, 7 
as in QMC. For comparison also results of high-temperature 
expansion (HTE) for J — Q are shown (dash-dotted line). 
The data for J — 0.4, 7 = 1 from Fig. M are not shown again. 



as a quasiparticle with a finite effective mass, exhibiting 
a cyclotron motion in B =/= 0. The latter behavior should 
lead to Xd(T — > 0) — > — 00 ||. Numerically it is easiest 
to test this conjecture for a single hole and 7 = (also 
true for J = 0). Namely, the QMC has no sign problem 
at 7 = 0, so that error bars are only due to the finite MC 
sampling. 

Results in Fig. || are quite interesting even for 7 = 
0. At J = a monotonous increase of \xd\ is observed, 
diverging as T — > M. It can be explained as a gradual 
transition from a hole in a random spin background to a 
well defined quasiparticle, i.e. the ferromagnetic polaron, 
at T — > 0. The situation is more complicated for J > 0. 
It is plausible that the difference to the J = case shows 
up at T < J, where the AFM short-range correlations 
appear. Relative to J = 0, a spin ordered state blocks 
the loop motion of holes, necessary for finite diamagnetic 
Xd- The effect is thus first a decrease of |xd(r) | with 
decreasing T, as seen in Fig. |9[ Turnover to the diverging 
diamagnetic Xd should happen only when the coherent 
quasiparticle is formed at Td <C J- Td should scale with 
the inverse of the quasiparticle effective mass I /to* . It 
is known that to* can become very large for the extreme 
7 = case, in particular for larger J. This explains why 
we cannot reach the coherent regime for J = 0.4 t even at 
T = 0.2t (Fig. ||), while the downturn is indeed observed 
for J/t = 0.1, 0.2, 7 = (Fig. |). 

At 7 > 0, the results are qualitatively different. The 
most pronounced effect is the change into a paramag- 
netic Xd for T < T p . The width of this T window is 
quite large. In fact within the FTLM and QMC data 
we are unable to locate the reentrance temperature Td 
into the diamagnetic response, although the latter is ex- 
pected B . An argument for the low value of Td can be 
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FIG. 10: Xd vs. Ch for several T and J = OAt. The last graph 
contains also a 4 th order HTE result (dotted). Canonical 
(dots) and grand-canonical (line) values for all ch are obtained 
within the FTLM with N = 16, where X12 is used to calculate 

Xd- 



given in terms of a very shallow energy minimum which 
defines the quasiparticle dispersion near the ground-state 
of a hole in an AFM within the t-J model [pT[ , hence the 
quasiparticle looses its character already at very low ex- 
citation energies. Still the paramagnetic response in the 
window Td < T < T p remains to be explained. 

Let us finally discuss also results_for \d for finite dop- 
ing Ch > 0, as presented in Fig. 
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to interpret is that of a nearly empty band, i.e. cy L > 0.7, 
where \<i is diamagnetic and nearly independent of T . In 
this regime the electron system is dilute and strong cor- 
relations are unimportant, hence Landau diamagnetism 
is expected. At moderate temperatures T > J and for 
an intermediate-doping regime, 0.2 < Ch < 0.7, Xd is 
dominated by a paramagnetic response with a peak at 
approximately Ch = 1/2. As consistent with results at 
low doping, there is a weak diamagnetism at Ch < 0.2 
and T > T v , while the paramagnetic regime extends 
to Ch = for T < T p . For low temperatures T <C J 
quite pronounced oscillations in Xd(ch) appear and can 
be partly attributed to finite-system effects. 

Certain aspects of the above results for Xd(ch) can be 
understood using the HTE. One is the asymmetry be- 
tween Ch — * and c/, — ► 1 at higher T > t. In the lowest 
order of the HTE, only hopping of the electrons around 
a basic plaquette loop has to be considered. The signal 
on the Ch — > side of Xd(ch) is thus reduced by a factor 
of 2 3 = 8 against the c/, — ► 1 side since in the first case 
only the plaquettes with ferromagnetically aligned spins 
contribute. 



algorithm solves some serious drawbacks of the MC meth- 
ods. It is indeed very efficient in cases where there is no 
sign problem, e.g. the anisotropic model 7 = at Nh = 1. 
Within the WL approach it is also very easy to formulate 
and measure certain responses like the orbital suscepti- 
bility Xd- Still the method suffers from a sign problem 
(for 7 > 0) even for a single hole Nh = 1 in an AFM 
(though not in a background of no holes Results 
are thus in practice limited to T > J/3 for the isotropic 
case 7 = 1. On the other hand the FTLM has no minus- 
sign problem but rather limitations due to small systems 
which can be studied. These are even more pronounced 
in cases with B > where the translational symmetry 
is lost. It is an interesting observation that within the 
FTLM, the limiting temperature Tf s for the t-J model is 
in most cases quite close to the lowest T reached by the 
QMC method. 

We have presented several results for thermodynamic 
quantities, i.e. energy E, specific heat C, entropy S and 
spin susceptibility Xs, as a function of T at low doping. 
At T > Tf s all results are consistent with the picture 
where holes introduced into the AFM behave as inde- 
pendent (nondegenerate) particles. The perturbation in- 
troduced into the AFM by holes is quite large even at 
lowest T < J, in particular visible from AC and AS, 
consistent with experiments in cuprates 63| . 

Results for the orbital susceptibility \di now obtained 
also for much larger systems using the QMC method, 
confirm the preliminary FTLM results || , indicating an 
anomalous paramagnetic response at low doping in an 
intermediate window of temperatures Td < T < T p ~ J 
(for the isotropic model 7 = 1). In fact within our nu- 
merical studies it is hard to reach the lower end of this 
window, meaning that Td < J/3. Still, the reentrance 
into a diamagnetic behavior is expected from theoretical 
arguments on the existence of a well defined quasiparti- 
cle at T — > 0, as well as from more reliable QMC results 
for the 7 = case ||. The paramagnetic response at 
intermediate T can be viewed also as an extension of a 
more pronounced Xd > regime observed at finite dop- 
ing 0.2 < Ch < 0.5 at all T. The explanation can thus 
go in the direction proposed by Laughlin || [27j], that 
at low doping Ch — > we are dealing with quasiparticles 
(with a diamagnetic response) , being a bound composite 
of charge (holon) and spin (spinon) elementary excita- 
tions. The binding appears however to be quite weak 
and thus easily destroyed by finite T or Ch, enabling the 
independent response of constituents, which apparently 
is paramagnetic. 



V. CONCLUSIONS 

Let us first compare both numerical methods used in 
the analysis of the t-J model at low doping. The QMC 
method allows the studies of larger systems and the loop 
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